#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include "dictionary.h"
#include "iniparser.h"

#include <gsl/gsl_math.h>
#include <gsl/gsl_errno.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_odeiv2.h>

#include "params.h"
#include "const.h"
#include "output.h"

#ifndef _REC_H_
#define _REC_H_

/* use log(a) instead of a or z as the independent varibale */
int ODE_rec( double lna, const double y[], double f[], void *params );

void Saha_Xe( cosmo_params *cp, double Tb, double *Xe );

void ODE_rec_test( cosmo_params *cp );

void do_rec_calculation( cosmo_params *cp );


#endif  // REC_H_
